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ABSTRACT 



We present 3D hydrodynamical simulations of a jet launched from the sec- 
ondary star of a binary system inside a proto-planetary nebula. The secondary 
star moves around the primary in a close eccentric orbit. From the gasdynamic 
simulations we compute synthetic [N ll] A 6583 emission maps. Different jet axis 
inclinations with respect to the orbital plane, as well as different orientations of 
the flow with respect to the observer are considered. For some parameter combi- 
nations, we obtain structures that show point- or mirror- symmetric morphologies 
depending on the orientation of the flow with respect to the observer. Further- 
more, our models can explain some of the emi ssion distributio n asym metries that 
are summarized in the classification given by ISoker &: Hadarl (120021 1 . 



Subject headings: ISM: jets and outflows — planetary nebulae: general — stars: 
AGB and post-AGB — methods: numerical — binaries: general 
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Introduction 



The shell of ionized gas coming from the ejected envelopes of the very late stage in the 
death of an intermediate initial mass (M,, < 8 M0) star is called Planetary Nebula (PN). 
Since the onset of high spatial resolution observational facilities, such as the HST, optical 
images of PNe have revealed a rich variety of morphologies, highly collimated structures 
and small scale features. These nebulae have been classified according to their large-scale 
morphologies as: bipolar, circular, elliptical or irregular. In particular, bipolar PNe are 
axially symmetric, having two lobes with an equatorial waist between them. 



Soker fc Hadarl (120021 ) proposed a classification of PNe depending on their departure 



from axisymmetry. The departures from axial symmetry can be point-symmetry or mirror- 
symmetry, with respect to the equatorial plane (which separates the two lobes). 

In the last decade, the point-symmetric shape of some PNe and proto-planetary neb- 
ulae (PPNe) has b een successfully explained in terms of collimated outflows. For example, 
Lee fc Sahail (120031 ) carried out axisymmetrical s i mulat ions for m odeling the morph ology and 
emission of the PPN CRL 618. I Velazquez et al.l (120041 ) (also see lRiera et al.ll2004j ). modeled 
the multiple shocked structures, emission and kinematics of the PPN Hen 3-1475 in terms 
of a precessing l e t with a periodic velocity variation, using 3D hydrodynamical simulations. 

elazquez et al.l (120071 ) model ed the thermal r adio-c ontinuum emission of the PN K3-35 as 
a continuous, precessing jet. iGuerrero et al.l (120081 ) present an observational and numeri- 
cal study of the PN IC 4634, showing that a variable velocity, precessing jet is the origin 
of its point-symmetric morphology. Some features observed in PNe (i.e. low excitation 
knots/cluraps moving at supersonic velocities) could be produced by cosmic "bullets" (e.g.. 



Raga et all 120071 : IPennis et al.l 120081 ). 



The presence of a companion to the AG B star is thought to be required in order to pro- 

he p resence of binary 



duce jets and point-symmetric morphologies (IBalick fc Frankll2002r). 



star system progenitors of PNe was initially proposed by lBond et a. 



1978 



, who discovered 



that UU Sge, the central star of PN Abell 63, is an eclipsing binary. iFabian fc Hansen] (119791 ) 
suggested that the particu l ar sha pe of the Helix nebula can also be due to the action of a 
binary system. iLivio et al.l (119791) sugges t ed tha t non-spherical PNe host close binaries. Fur- 
thermore, the survey of iDe Marco et al.l (120081 ) shows strong evidence that at least 10-15% 
of PNe are composed of binary systems with short orbital periods (< 3 days). Eventhough 
this fraction is not high enough to be consistent with the large fraction of PNe that appar- 
ently have been shaped by binary interactions, this survey is not able to detect large orbital 
periods and therefore the binary fraction of 10-15% has been considered a lower limit. 



Some mechanism invoked to explain the shaping of bipolar PNe involve the action of jets 
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launched by the primary or secondary star (e.g.. iMorrisI 119871 . 119901 : iSokerl 119921 : iLiviol 11993 
Soker fc Rappaportlboooh . In several models, the jets are supposed to form when a main 
sequence star (or a white dwarf) accrete material from the AGB (or post-AGB) star, forms 



an ac cretion disk and blows two jets (e.g.. lFrank fc Blackmanll2004l : lFrankll2006l : iFrank et al. 



20071 ). Hydrodynamical simulations have shown that jets (blo wn by the AGB star or by a 
compa nion) can account for some of the observed morphologies. iGarcfa-Arredondo fc Frank 
( 12004| ) carried out 3D numerical simulations of the interaction of the stellar winds from the 
components of a binary system. They found that bipolar morphologies as narrow-waisted 
nebulae can be obtained after the interaction of a slow AGB wind, from the primary star, 
with a coUimated fast stellar wind produced by the secondary. 

In this work, we explore the influence of having a binary system in order to reproduce 
the morphology and brightness distribution of bipolar PPNe. We investigate (3D hydrody- 
namical simulations) if both point- and mirror- symmetric PPN could be explained in terms 
of an interacting jet (source in orbital motion) with an AGB wind. Also, we explore if the 
observed morphology of the PPN is affected by the orientation of the flow with respect to 
the observer. The jet is ejected by the secondary star (a star in the main sequence), which 
is on orbital motion around the remnant of the AGB star. For simplicity, the AGB wind is 
considered isotropic, i.e., it has a density distribution given by oc and a constant wind 
ve locity (see section 3) . We consider a clos e eccen tric binary system, calling the suggestions 
of ISoker et aP (|l998h . ISoker fc RappaportI tood ) and lSoker fc RappaportI (|200lh . 



2. Modeling the orbital motion of a binary system 



We investigate the interaction of a jet with the AGB slow wind in a binary system 
scenario, where the primary is the remnant AGB star. Observa t ions show that many PNe 
seem to have binary sources fjvan Winckell l2003l : IZijlstral 120071 : iDe Marco et al.l |2008| ) , in 
which the two stars can be in different evolutionary stages according to their initial masses. 

Both stars orbit around the center of mass in stable elliptical paths. The Keplerian 
velocity of the reduced mass is corrected by the factor Mp/ {Mg + Mp), in order to obtain the 
orbital velocity of the secondary star. The orbital plane (the xy- plane) is such that both foci 
of the ellipse are located along the x-axis, i.e. both the periastron and apoastron lie along 
this axis. With this configuration, the orbital velocity achieves its maximum magnitude in 
the periastron, toward the positive y-axis (see panel (a) of Figured]). 



van Winckell (120031 ) gives values for the semi-major axis of the order of 1 AU, while 
PNe usually have physical sizes of ~ 0.1 pc (~ 10^ AU). As we are interested in modeling 
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the large scale structure of PNe, the spatial resolution of our numerical simulations is not 
high enough to resolve the scale of the orbital motion. However, our numerical simulations 
do resolve the temporal scale of the Keplerian motion, so that we can include the orbital 
speed as an additional component of the outflow velocity. In panels (b) and (c) of Figuredl 
we display the actual scheme used, i.e. a fixed jet, centered in the computational domain, 
where the orbital speed is added to the jet velocity. Two snapshots are shown for the cases 
where the orbital speed is maximum -at the periastron, panel (b)- and minimum -at the 
apoastron, panel (c)-. 

A detaile d exploration of the pa rameter space (semi-major axis, eccentricity, etc) is 



carried out by lHaro-Corzo et al.l fl2009l . in preparation). In this work, we only show the most 
relevant results, choosing Mp + Ms = 1.25 M0(with Mp = IMq) and a semi-major axis of 
3.2 AU (i.e., an orbital period of 5 years). 



3. Initial setup for the numerical simulations 



The 3D numerical simulations were carried out with the yguazu-A code (IRaga et al. 



2000l ). The yguazu-A code integrates the gas-dynamical equations together with a system of 
rate equations for the atomic/ionic species: HI, HII, Hel, Hell, Helll, CII, CHI, CIV, NI, Nil, 
NIII, 01, OH, OIII, OIV, SII and SHI. With these rate equations, a non-equilibr ium cooling 



functi on is computed. The reaction and cooling rates are described in detail by iRaga et al. 



(120021 ). The gas-dynamical equations are integrated with a second ord er accurate tech nique 



(in time and space) employing the "flux- vector splitting" algorithm of Ivan Leerl (119821 ). 



Our numerical simulations use a 5-level, binary, adaptive Cartesian grid, in a computa- 
tional domain of 256 x 256 x 512 pixels (at the highest grid resolution) along the x-,y-, and 
z- axes. The size of the computational domain is of 2 x lO^^cm along the x- and y-directions, 
and 4 x 10^^ cm along the z-direction. 

A jet is injected into a circumstellar medium, which has been swept up previously by a 
dense and slow wind from the AGB central star. Even though the circumstellar medium in 
bipolar PNe is aspherical (with higher densities in the equator than along the bipolar axis), 
for simplicity and in order to only study the influence of the binary orbital motion, we 
have "turned off" the AGB influence imposing on all computational domain the analytical 
solution of an isotropic AGB wind. Then, the wind density distribution is given by: 



w 



4 TT v„ 
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where is the mass loss rate of the AGB wind, Vw is the AGB wind terminal velocity, 
and r is the distance from the primary star. We assume that the AGB wind is initially 



neutral, with = 2 x 10"^ yr" 



1.5 X 10^ cms ^ and a temperature = 100 K. 



A decreasing temperature vs. radius dependence (resulting from the adiabatic expansion of 
the spherically diverging wind) is generated as the time-integration progresses. 

The jet is injected at the center of the computational domain within a cylindrical volume 
with radius and length Ij with values = = 4 x 10^^ cm. (which are equivalent to 
5 pixels in the highest resolution grid). The jet temperature was set to 10^ K. The jet 
velocity is = 1.7 x lO^cms"^ and the mean atom + ion number density is constant with 
rij = 5 X lO^cm"^. These values correspond to a Mj = 1.5 x 10~^ M0 yr~^ mass injection, 
which could seem high. Due to the limited numerical resolution we are forced to impose the 
jet on a region larger than its true size. H owever, the total injec ted mass after 400 yr is 
1.2 X 10~^ M0 , which looks like reasonable (j Velazquez et al.l 120071 ). 



With the par ameters listed above, the je t is i n the domain of "strong jet" according to 
the description of ISoker fc RappaportI (I2OOOI ) and iGarcia-Arredondo fc FrankI (|2004J ). with 
a parameter x = M^v^j/MjVj = 0.2. A similar dynamical behaviour could be achieved 
taking into account a less dense jet but with higher velocity, keeping constant the jet ram 
pressure pjVj and the parameter x- However, we prefer to employ a dense jet with a low 
velocity, being the last one comparable with the maximum orbital speed, in order to favor 
the contribution of orbital motion to the global morphology of the object. Based in these 
initial conditions, we consider six runs labeled as Rai and Rbi, where i=l corresponds to the 
case in which the jet axis is perpendicular to the orbital plane, i=2 to a tilt of 20° between 
the jet axis and the orbital axis, on the x^-plane, and i=3 to a tilt of 20° on the yz-plane 
(the orbital axis is parallel to the 2-axis, and the x-axis contains the foci of the elliptical 
orbit, see §2). Runs Rai have an eccentricity e = 0.73, and runs Rbi have e = 0.94. The 
parameters of the runs are summarized in Table [TJ 

Finally, from the temperature and atomic/ionic/electronic number density distributions 
computed from the numerical time-integration, we calculate the emission line coefficients of a 
set of permitted and forbidden emission lines. The forbidd en line [N i l] A 65 83 was computed 
by solving 5- level atom problems, using the parameters of iMendozal (119831 ). 



4. Results 



We find that many of our jet/AGB wind interaction models have morphologies that 
depend on the projection, and that both point- and mirror- symmetric PNe can be obtained 
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from the same model (depending on the position of the observer). This is illustrated in Figure 
|2l where we plot the resulting synthetic [N ll] A 6583 emission maps at an integration time of 
400 yr. The group of six panels correspond to the runs Rai. The yz- and xz-projections of 
run Ral show mirror-symmetric morphology. On the other hand, while the yz- projection 
of run Ra2 has a mirror-symmetric morphology, the a;2;-projection of this model has a point- 
symmetric morphology. A similar switch between point- and mirror-symmetric morphologies 
is seen in run Ra3. Also, we analyze the asymmetries of the synthetic [Nil] A 6583 emission 
maps for runs Rai produced by the inclination (40°) of the orbital plane (x2;-plane) with 
respect to the plane of the sky. The intensity maps preserve the mirror- and point-symmetric 
morphologies discussed above (for the case in which the 2;-axis lies on the plane of the sky). 
Furthermore, we find that the morphologies of the maps obtained from runs Rbi are similar 
to the maps from Rai, although in general, we note the presence of small-scale internal 
structures which are brighter for runs Rbi than runs Rai. Those are because of the larger 
eccentricity of the Rbi runs. The velocity variability of the jet (resulting from the orbital 
motion), produces these internal substructures in both lobes of the Rai and Rbi maps. 
In order to see the influence of the increasing eccentricity, in the figure[3] we display the 
substraction of the synthetic [N ll] A 6583 emission maps corresponding to Rbi- Rai runs. 

Finally, we have obtained the [N ll] A 6583 luminosities for the top and bottom lobes for 
all of the computed maps. Table [2] gives the [N ll] A 6583 luminosities for both lobes and 
all runs. In general, we note that runs Rbi have larger [N ll] A 6583 luminosities than runs 
Rbi (of the order of a 10%), as a result of the more eccentric orbital motion of runs Rbi. 
The exception is the run Rb3 because the total luminosity measured on the top lobe is the 
same of run Ra3, while the bottom lobe has less total luminosity. This is due to the gas 
material of top lobe has larger velocities than the bottom lobe gas. In the top lobe, the 
orbital motion speed is added to the jet velocity, while the opposite occurs for the bottom 
lobe. For the cases where the resulting shape has mirror-symmetry, the luminosities from the 
right and left sides of the lobes have also been computed. Comparison of total luminosities 
between both (top and bottom) lobes show a similar value for runs Ral, Ra2, Rbi and 
Rb2. However, analyzing the xz- and yz- projection cases, an asymmetry in the brightness 
distribution is observed for the ?/2;-projection, in general, with the right side being more 
luminous, consequence of the orbital motion. In runs Ra3 and Rb3 we obtain appreciable 
[Nil] A 6583 luminosity asymmetries between the top and bottom lobes (see figure |3]). 

Point- and mirror-symmetric morphologies are simultaneously obtained if the jet axis 
has a certain inclination with respect to the orbital plane. Interestingly, the classification of 
a PN belonging to one or other group seems to depend on the point of view of the observer 
or on the projection. There are significant differences in the total luminosities of both lobes 
between runs Ra2 and Ra3 and similar behavior is observed for runs Rb2 and Rb3. Run 
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Rb2 produces top and bottom lobes with similar total luminosities, while run Rb3 generates 
a luminosity difference of at least 36% between the top and bottom lobes. This is due to the 
fact that the jet axis is tilted 20° on the yz-plane in the case of run Ra3. In this case, the jet 
inclination contributes to the elliptical motion, being its effect larger for the top lobe when 
the secondary star passes by the periastron, where the orbital motion velocity is maximum. 
The total jet velocity (the intrinsic constant jet velocity ~ 1.7 x lO^cms"^ plus the eccentric 
orbital motion velocity) is also maximum at the periastron, producing strong shocks and 
increasing the emission for the top lobe (the opposite occurs for the bottom lobe). 



5. Discussion and Conclusions 

In this work, we carried out 3D hydrodynamical simulations with the yguazu-A code 
in order to explore the binarity mechanism for shaping bipolar, point- and mirror-symmetric 
PN. We consider a collimated outflow (jet) launched by the secondary star of a binary system, 
located at the center of a PN. The primary star is the source of the AGB wind, in which the 
jet is embedded. The two stars have a close, very eccentric orbit, resulting in a maximum 
orbital velocity of the order of the jet velocity. Because of this, the orbital velocity has a 
strong effect on the jet dynamics, and strongly modifies the emission line maps predicted 
from the computed flows. 

We have generated [N ll] A 6583 emission maps considering different inclinations of the 
jet axis with respect to the orbital plane and also different tilts between the orbital axis 
and the line of sight. We find that the predicted maps of some runs exhibit (see figure 
[2]) either point- or mirror-symmetric substructures depending on the orientation of the flow 
with respect to the observer. Also, the maps exhibit substructures with different luminosities 
for the two lobes and with side-to-side asymmetries within the same lobe. 



In this sense, ISoker fc Hadarl (120021 ) built a classification of PNe based on their departure 
from axisymmetry considering both morphological and brightness distribution criteria. From 
Figure[2]and tabled we n ote that runs RaS and Rb3 would correspond to class 2 (i.e. unequal 



intensity of both sides) of ISoker &: Hadarl (120021 ) . because top lobes are brighter than bottom 



ones. Into this class could be also included maps of runs Ral and Ra2, for example, because 
in their |/2;-projection display an asymmetry in the right and left brightness distribution. The 
^/^-projections of Ral and Rbl could be also included in the class 4 (i.e. the jets are bent 
to the same side, therefore there is a mirror-symmetry), because the small-scale structure 
observed on the right side of both lobes. 



From our study we conclude that a jet ejected from a companion to the AGB star 
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(which gives rise to a PN) in a highly eccentric orbit leads to morphologies with either point- 
or mirror-symmetries (in some cases both, depending on the inclination to the observer). 
Also, the bipolar nebular structure resulting from the jet/AGB stellar wind interaction has 
brightness asyrametri es with contrasts similar to the asymmetries observed in some PNe 
( ISoker fc Hadarl l2002l ). We must note that, except for the eccentricity, the orbital pa- 
rameters we used are commonly found in last evolution state of low m ass and small post 



AGB stars, which can be a close binaries without a common envelope (Ivan Winckell 12003 



Frankowskil l2004l : IVan Winckel et al.ll2009l : Ide Marcol 120091 ) . The eccentricities we used are 
uncommon for post- AGB stars, h owever, there are reports of highly eccen t ric orbits in other 
binary systems, such as Ba stars ( van Winckell l2003l : IJorissen et al.lll998l ). Ide Marcol ( 20091 ) 
argued that for post-AGB stars that have been monitored spectroscopically, low eccentrici- 
ties have been found for systems with intermediate periods (between 100 to 1500 days), but 
for periods > 1 yr the eccentricities can be close to unity). A larger sample of post-AGB in 
binary systems is needed to assess the feasibility of our model. 

It is interesting that we can generate models which produce either point- or mirror- 
symmetric emission line maps, depending only on the orientation of the flow with respect 
to the observer. For other choices of flow parameters, however, it is easy to obtain models 
which produce point-symmetric structures regardless of the orientation of the flow with 
respect to the observer (e. g., b y having a precessing jet from a source with low orbital 
velocity, iMasciadri fc Ragall2002l ). Therefore, while a large range of parameters will produce 
point-symmetric nebulae, we find that it is necessary to have both an appropriate parameter 
combination and an appropriate orientation of the flow (with respect to the observer) in 
order to produce mirror-symmetric structures. This result might be consistent with the fact 
that many more bipolar PNe are observed to have point-symmetric structures (as opposed 
to mirror-symmetric structures). In order to decide whether this is indeed the case, a more 
detailed study limiting the parameter space relevant for nebulae resulting from jet/AGB 
wind interactions will have to be made, as well as a study of the statistics of the resulting 
flow morphologies observed at random orientations. 
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Fig. 1. — (a) Scheme of a vertical jet in orbital motion. The center of mass is located in 
one focus of an elliptical orbit and the secondary star orbits in an elliptical path around the 
center of mass (which lies close to the primary star). The foci are along the x-axis. The 
xy-plane is the orbital plane, (b) and (c) are the true scheme employed in our simulations, 
where the jet is located at the origin. The orbital speed is added as an additional component 
of the jet outflow. 
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Ral Ra2 Ra3 




Fig. 2. — Synthetic [Nil] A 6583 emission (erg s~^ cm~^ sr~^) at 400 yr for runs Ral, Ra2 
and Ra3. The emission is on log gray-scale. Axes are in units of 10^'' cm. [See electronic 
edition for a color version.] 
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Table 1 


: Run paramet- 


ers. 


run 


ecc 


angle (°) 


plane 


Ral 


0.73 







Ra2 


0.73 


20 


xz 


Ra3 


0.73 


20 


yz 


Rbl 


0.94 







Rb2 


0.94 


20 


xz 


Rb3 


0.94 


20 





Table 2: Luminosity of top and bottom lobes L/Lq(10 ^). Each lobe has been parted as left 
and right side. We ignore the central emission (20 pixels from each lobe). 

Run Bottom xz yz Top xz yz 

total left right left right total left right left right 

Ral 3.4 1.7 1.7 1.3 2.1 3.4 1.7 1.7 1.3 2.1 

Ra2 6.9 - - 3.2 3.7 6.9 - - 3.2 3.7 

Ra3 5.6 2.9 2.7 - - 6.9 3.6 3.3 - 



Rbl 
Rb2 
Rb3 



3.8 
7.6 
5. 



1.9 



2.6 



1.9 



2.4 



1.5 
3.4 



2.3 
4.2 



3.7 
7.4 
6.8 



1.^ 



3.6 



1.9 



3.2 



1.5 
3.3 



2.2 
4.1 
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Fig. 3. — Synthetic [Nil] A 6583 emission at 400 yr. for the substraction Rbi- Rai. Same 
scale as figured [See electronic edition for a color version.] 



